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Abstract. It is well known that a large fraction of galaxies have cuspy luminosity 
profiles in their central regions, at least within the observational resolution. In such 
cases, the often used, simplified, local approximation for the dynamical friction brak- 
ing classical term fails when the massive satellite moves through the inner parts of the 
galaxy, although the scattering integral still converges for phase space distribution sin- 
gularities that are not too sharp. Here we present preliminary results of our work aiming 
at finding better and more reliable results from the integration of motion of massive ob- 
jects (globular clusters) in galaxies where the density diverges to the center in a power 
law form, with exponent greater than -2. 



1. Dynamical friction 

The frictional deceleration of the motion of a test particle of mass M placed, at time 
t, at r is given by the cumulative effect of a background of target stars of mass m at 
varying the impact vector b, which can be expressed as 



d\ M 2m r r v M -v m |v M -v m | 3 3 

—T = -T7^\ /(r + b,v m ) d^d'b, (1) 

where f(r + b, v OT ) is the stellar background distribution function (DF) at position r' = 
r + b. Even for a DF inspace and velocity spherical simmetry, i.e. f(r', v m ), the integral 
in Eq.(Q]) is very complicated whenever the test particle is out of the origin (r + 0). 

To overcome, partially, this problem, the common approach consists in taking out 
of the integral in Eq. (Q]) the DF, assumed as independent of b and isotropical in the 
velocity space, thus yieldin g to what is often referred as local approximation formula 
Binnev & Tremaind Il987h 
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-4nG 2 Mp(r)F(r,v M )lnA^-, (2) 



the function < F(r, vm) < 1 being the local fraction of target stars slower than the test 
object (v m < vm). In the local approximation, the mass density, p, of field stars in Eq. 
(0 is assumed at its local value, i.e. that evaluated at the position (r) of the test particle; 
In A is the usual Coulomb's logarithm. 

Now, whenever the test particle is significantly off center with respect to the stel- 
lar system (star cluster, galaxy, etc.) wherein it moves, the local expression gives (for 
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Figure 1 . The time evolution of the distance to the galactic center of a massive 
test object in the case of initially circular orbits (regular curves) and radial orbits 
(oscillating curves) of same initial orbital energy, for two values of y, as labeled, and 
two values of the object mass (upper curves refer to M — 10 6 M . lower curves to 
M - Ax 10 6 Mq). 



symmetry reasons) a quite acceptable approximation; on the contrary, it loses its valid- 
ity in the neighbourhood of the centre. In this case the local approximation is clearly 
an overestimate of the actual dynamical friction, because it corresponds to weighting 
the contribution of the gravitational encounters with objects at a given distance from 
the test particle not with the, correct, density of target stars at that distance but, rather, 
with the density of targets evaluated at the location of the test particle itself, that is 
maximum at the origin of any self-gravitating system. This overestimate is particularly 
relevant when dealing with cuspy galaxies, where the spatial density of stars diverges 

at the galactic center, making the local approximation unviable. 

To solve this probl em, we follow the approach presented in IVicari et all (120071) and 
deeply discussed in lCapuzzo-Dolcetta & Arca-Seddal ((201 lh ). which is based on taking 
into account that at the origin (r = 0) and in the hypothesis of isotropic DF the integral 
in Eq. £T|) simplifies into the form 



(d\M\ 87rm C bmax f y M -y„, , 

I dt I cen M + lfl J hmm J Ym ! + b X">r M l 

which can be performed numerically, although with the care due to the double singu- 
larity in phase space. Once that the integral in Eq. © is evaluated, an estimate of the 
dynamical friction deceleration along the test object trajectory is given by an interpola- 
tion between the local approximation as given by Eq. [2] and the central value given by 
Eq. ©, in the form 

d\M , Jd\ M \ r , . .Jd\M\ ... 

^ = ^hrL +v - mi hrL- (4) 

where the interpolation function < p(f) < 1 is such that p(0) = 1. We found 
p{r) = e~ r l R (where R is a proper length scale) as a good interpolation formula. We 
have applied this to the orbital evolution of massive objects in a spherical potential 
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resembling that of the Milky Way (disk excluded), as given by the s um of three c ompo- 
nents: i) a cuspy inner bulge, represented as a Dehnen' s y-model (|Pehnenlll993 ). ii) an 
extended bulge-halo represented as a Plummer sphere (Plummet! 1911) and, iii) a dark 
matter halo in the Navarro, Frenk and White form dNavarro et al. 1996 ). 



2. Some results 

The comparison between the efficiency of dynamical friction on circular and radial or- 
bits of same energy is particularly important because they constitute (for any given 
orbital energy, in a spherical potential) the lower and upper limits of the expected dy- 
namical braking. In Fig. [T]we show, for two values of y in the Dehnen's model, the 
time evolution of the distance to the center for circular orbits starting from r = 1 kpc 
from the galactic center and for the radial orbits of same energy, for 2 values of the test 
object mass, M = 10 6 M G and M = 4 x 10 6 M . The ratio of the radial to circular 
dynamical friction decay times, i.e. the time needed to lose almost all the initial orbital 
energy, is about 0.04 for both y = 1/2 and y = 1, indicating how efficient is dynamical 
frictional braking through galactic, cuspy, centers and, so, how important is its careful 
treatment. 
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